Resonantly Tunable Majorana Polariton in a Microwave Cavity 
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We study the spectrum of a one-dimensional Kitaev chain placed in a microwave cavity. In the off-resonant 
regime, the frequency shift of the cavity can be used to identify the topological phase transition of the cou- 
pled system. In the resonant regime, the topology of the system can be controlled via the microwave cavity 
occupation and, moreover, for a large number of photons (classical limit), the physics becomes similar to that 
of periodically-driven systems (Floquet insulators). We also analyze numerically a finite chain and show the 
existence of a degenerate subspace in the presence of the cavity that can be interpreted as a Majorana polariton. 
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Introduction. — Majorana fermions, or half-fermions, have 
recently attracted tremendous attention as building blocks of 
a topological quantum computer Uj-]4J. They are predicted to 
emerge as excitations in several solid-state systems, such as 
genuine two-dimensional (2D) p-wave superconductors Q, 
or induced via superconducting proximity effects in topolog- 
ical insulators S |6] |71 or one-dimensional (ID) wires ||8l- 
[T2| as boundary zero modes. They obey non-Abelian statis- 
tics, both in 2D L13J and ID [14J, allowing for implemen- 
tation of certain (nonuniversal) gate operations required in 
quantum-computational schemes via braiding of the Majorana 
fermions. Moreover, due to their highly nonlocal character, 
the qubits built out of Majorana fermions are insensitive to 
local parity-conserving perturbations, thus making them po- 
tentially robust against physical noise errors 121 [TSl [T6l . 

The existence of Majorana states, however, has been shown 
not only for static systems, but also for driven systems, as 
zero modes in the so-called Floquet Hamiltonians lfT7l4T9ll . or 
even, more recently, in dissipative systems EOl . These pro- 
posals are based on the idea that driven systems can have 
different topology from their static parents, albeit these are 
not actual ground states. The ground-state perspective can, 
however, be used even for the driven models, if instead of a 
classical driving one considers the driving as having its own 
quantum dynamics, such as in a high-Q electromagnetic cav- 
ity. In particular, ID microwave cavities have been proven ex- 
tremely successful in reaching the so-called strong-coupling 
regime between photons and different types of qubits |I 2T1|291 . 
with very large Q factors and a high degree of control. More- 
over, strong coupling between atomic gases (e.g., in a Bose- 
Einstein condensate) and optical cavities have been achieved 
with a high degree of coherence BOll . In this paper, we com- 
bine cavity QED with Majorana physics in a ID lattice model. 
We analyze the spectrum of a p-wave lattice superconductor 
(Kitaev chain) coupled to a microwave cavity and examine the 
topology of the combined system. The spectrum is studied us- 
ing Dicke-Iike Hamiltonian, in both the off- and on-resonant 
regimes. 

The model. — In Fig. [T| we show a schematic of the sys- 
tem under consideration: a Kitaev chain inserted in a ID mi- 
crowave cavity. The physical sites are depicted by black ovals, 
while the constituent Majoranas are shown inside by solid red 




FIG. 1. A sketch of the system. In blue is the superconducting mi- 
crowave cavity, which supports quantized electromagnetic modes. 
The Majorana chain yj 2 • ■ • jlf is depicted inside the cavity with the 
relevant parameters renormalized due to the coupling to the cavity: 
Chemical potential ji ^ p + 6fi and hopping amplitude f — » t + 6t 
(with the hats denoting fluctuating quantum contributions). 



circles. The hopping parameter t and the chemical potential fi 
can be modified by the (quantum) cavity electric field E(x), 
engendering a physical coupling between electrons and pho- 
tons. The associated wavelength is assumed to be much larger 
than the ID chain, so that the electric field has no spatial de- 
pendence along the latter. 

The total Hamiltonian of the system, 'K - Hiu+ Him + H^h, 
consists of 



N-l 



HiD^-J] {tcl^cj + AcJ^jcJ) + H.c. - ^ c]cj , 

' N N '• 

Hint = y6 X ^1+1 ^-/^ + H-c. + a ^ c]cj (a"'" + a) , ( 1 ) 
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and i/ph = 0) era, where t is the hopping parameter, A is the 
p-wave pairing potential, jj is the chemical potential, a is the 
electron-photon coupling that shifts the chemical potential, 
and j6 accounts for changes in the tunneling Hamiltonian in 
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the presence of photons. While the hopping change in elec- 
tric field appears naturally, the modification of the chemical 
potential is unambiguous only in the presence of a supercon- 
ducting reservoir and is associated with supercurrents flowing 
from and into the nanowire on a time scale determined by the 
strength of the associated Josephson coupling. We are thus 
tacitly assuming that such Josephson oscillations are much 
faster than the relevant dynamics in the cavity, so that the ef- 
fective pairing A is essentially anchored at a constant value. 
For simplicity, we also assume that A is unaffected by the 
photon field, a (a') stands for the photon creation (annihi- 
lation) operator and u is the frequency of the corresponding 
photon mode (setting h - \ throughout). The electric field in- 
side the cavity (correspondng to its fundamental harmonic) is 
E{x) - By yjiolLcP-csin (nzlL)(a^ + a) |23|, with L being the 
cavity length, d is the distance between the center conductor 
and the ground conductors (see Fig. [TJ, and c is the capaci- 
tance per unit length. Note that the maximum of the electric 
field is at z = L/2, where we suppose to position the wire. 

Isolated wire. — We begin by recapping some of the known 
results on the Majorana physics in ID p-wa\e superconduc- 
tors. To keep the ideas transparent (but without loss of gen- 
erality), we assume in this section real f = A > and 
fi - 0. Performing the substitutions yj = cj + cj and 

- i(cj - cj), in terms of Majorana operators, (yp' - y'j, 
with / - 1,2, we can write the ID wire Hamiltonian sim- 



ply as Hid = -''2^=i' T^rj+p with y\ and dropping out 



entirely. This guarantees a degenerate ground state, which 
can be viewed as a qubit indexed by cr, = (c|,Cf - 1/2). 
Cf = (y[ + iyl/)/2 here is the nonlocal fermionic operator 
defining the corresponding zero-mode quasiparticle. 

The existence of the Majorana end modes can be identi- 
fied using bulk properties only, as they are a consequence of 
the topology of the Brillouin zone. Using periodic boundary 
conditions allows us to switch to the reciprocal space using 
Cj - Tjk Ck exp {-ikj)l ^In to write *Hid = ^^BdG(^), where 



HBdc(k) - -(t cos k + iu)tI - tsink 



(2) 



is the Bogoliubov-De Gennes Hamiltonian (restoring a fi- 
nite fi), and the pseudospin Tj. - (t^, rp, expressed 
in terms of the Pauli matrices, acts on the particle-hole 
basis (Q,ci^). One can diagonalize this Hamiltonian 
HBdG(^2 " UHk)HsdG{k)U{k) by Uj k) ^ exp(/gtTp, w hich 
gives //fidG = ^kTl, with ek — ^f-+iJ?-+ It/icosk and 
Ok = arctan [f sin^/(f cos^ -H ju)]/2. The topological invariant 
(winding number) that measures the number of edge modes at 
each end of the wire reads [,3 1 il : 



(3) 



and gives P - 0(1) for t < (t > \fi\), implying the absence 
(existence) of Majorana end modes for a system with open 
boundaries. We will now analyze this quantity in the presence 
of the cavity. 




FIG. 2. The dimensionless functions r{x) (inset) and r2(x) as a func- 
tion of X = \fj.\/t. There is a strong signature of the transition point 
X = 1 (i.e., \fi\ = t) in which becomes singular. 



Cavity-coupled wire. — The combined system Hamiltonian 
can be mapped to the well-known Dicke model, namely a set 
of spins interacting with the same photonic mode. We start 
with the bulk system, later focusing on a finite wire. In the 
continuum limit, we can write the interaction Hamiltonian as 
^^int - ctYik^kifl^ + ci)Tj^, with 6k - I + (J3/a) COS k, which re- 
sembles the interaction Hamiltonian between a photonic mode 
and a collection of spins: the Dicke model. To analyze the 
effect of this term on the spectrum, we consider both the high- 
frequency (w » f H- \fi\, i.e., much larger than the quasiparticle 
band width) and the resonant (|f - \fi\\ < cj < t + \fi\) regimes. 

Off-resonant coupling. — For a large cavity detuning, lj » 
t + |yu|, all the transitions are oflF-diagonal, and we can employ 
second-order perturbation theory using the Schriefifer-Wolff 
formalism. This means diagonalizing the Hamiltonian in a 
decoupled basis of the wire-cavity system via a unitary trans- 
formation, •K = exp (5)'?/ exp (-5) = H- [5 , TY] H- 0(S\ 
with S ' - -S being an anti-Hermitian operator. We leave the 
details for the Supplementary Material (SM) and only mention 
that, although in this situation the photon field cannot easily 
change the topology of the system, one can utilize the cou- 
pling to measure the topological phase of the system[32 | as 
well as the transition point. This is encoded in the cavity fre- 
quency shift 6(D - ANtiajcSf- Tin It) with 

r(/^/f) ^Y^Yi^l ^<^^'> * ^ J "^^^l ^<^^> ' (4) 

where (...) means averaging over the (mean-field) ground 
state of the ID system. The perturbative treatment is justified 
for finite-size chains with a y/N <«c co. Moreover, we note that 
besides the frequency shift, the cavity leads to a shift of the 
effective chemical potential, which is small in the perturbative 
regime (see SM). In Fig. [2] we plot r(x) and r2(x) = d^T/dx^ 
as a function of the chemical potential fi. We see that while 
r(x) varies smoothly, r2(x) becomes singular at the topolog- 
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ical transition point, the divergence being logarithmic in na- 
ture, i.e., r2(l + e) oc log|e|. The corresponding detection 
method for the transition would thus provide an optical al- 
ternative to the more conventional transport-based proposals 
such asthatofRef. f33ll. 

In order to experimentally resolve the frequency shift, the 
condition dcDjci) » Qr^ must be met, with Q being the qual- 
ity factor of the cavity. If we assume microwave cavities with 
Q ~ 10^ 1^, and taking t/u ~ 10"' and ay/N/oj ~ 10"^ 
we get d(L)/a) ~ 10"^ » Q"'. Note that even if the coupling 
parameter a to the individual bonds turns out to be small, by 
assuming a large enough chain (i.e., large A^) one can drasti- 
cally enhance the effective coupling a y/N to reach the strong 
coupling regime. Moreover, since the cavity is highly de- 
tuned from any transitions in the wire, no additional relaxation 
mechanism are introduced by this coupling. 

As an extension of our original model, Eq. ([TJ, we may con- 
sider an alternative interaction Hamiltonian, Him - g 
a' r^), where = (r^ ± /rp/2, which does directly affect the 
topology of the system. Note that we do not suppose here a 
rotating-wave approximation (see below). To the second order 
in the coupling strength g, we find: 



a a + 



1 



ki=k' 



(5) 



The first term here accounts for a cavity-induced shift of the 
chemical potential, while the flip-flop coupling in the second 
term leads to only a small mean-field chemical-potential shift 
in the ground state, similarly to the previous section (see SM). 
With this effective interaction, it is possible to switch between 
topologically trivial and nontrivial states by changing the pho- 
ton number (which is conserved). For example, neglecting 
the second term in Eq. (jSj), the first term leads to a shift in 
the chemical potential, n — > jt/eff - - (g^/aj){nph + 1/2), 
with = (a^a) (assuming henceforth that fj. > 0). If the 
system is in the nontopological phase in the absence of the 
cavity coupling (i.e., /j > t), for Wph > - (fi - t)(tL)/g ) 
the effective chemical potential satisfies yUgff < f, thus the 
system crosses over to the topological phase (we neglected 
the 1 /2 part since typically this transition would happen for 
Wph s> 1). This shows the attractiveness of such a modified 
Hamiltonian, which is a quantized analogue of the one pro- 
posed in Ref. |'34l|. This type of Hamiltonian can, however, 
be physically implemented only if the coupling to the cavity 
leads to changes of the pairing field. 

Resonant coupling. — In the resonant regime, it is possi- 
ble to change the system topology at the single-particle level 
of Eq. ([T]) depending on the cavity state. To see this, we 
first simplify Hamiltonian ([T]i by neglecting the "counter- 
rotating" terms that are off-resonant. Leaving the details of 
the derivation for the SM, we show here only the final result 
[in the rotated particle-hole basis as described below Eq. Q]: 
A^int = Yjk ^'^('''^'^ + <2"'"t^ ), where a-^ = aS^ sin (20^). This 
closely resembles the original Dicke Hamiltonian |35 1. How- 
ever, note that both the single-pseudospin splittings as well 




FIG. 3. Left: single-particle spectrum £(k) for zero photons in the 
cavity (or in the absence of the cavity). Right: single-particle spec- 
trum in the presence of the cavity coupling, E„(k), for difFererent 
number n of photons in the cavity. For the plots we used: fi = 0.3, 
a = 0.1, jS = 0, and ai = 0.9, all in units of t. Note that these 
parameters put us in the resonant regime. 



as the couplings are k dependent (as opposed to the model 
in the previous section). For each individual k, the Hamilto- 
nian reduces to the Jaynes-Cummings Hamiltonian, which is 
block-diagonal, with each block being a 2 x 2 matrix. Each of 
these blocks is associated with a given value of the conserved 
quantity Ck = + a'a, so that the 2x2 block acts in the 
subspace {Itt) ® |n),lijt) ® \n + 1)). This results in an eff'ec- 
tive two-band Hamiltonian, which we can write in terms of a 
pseudospin s = (sx, Sy, s,) (Pauli matrices) as follows: 



1 y/nal 



(6) 



giving the spectrum E„(k) - ± ^(lo - e^)^ + n(a^y-/2. Note 
that we neglected a constant shift AE„ = (n H- 1 /2)w of the 
two bands for a given - n. For each block labelled by 
n, the spectrum has a band gap given by the strength of the 
pseudospin-photon interaction times This dramati- 
cally alters the initial spectrum e^, as n increases, which is 
depicted in Fig. |3] 

The single-particle picture described above is intimately 
related to the recent developments on the topology in 
periodically-driven systems |[T7Hl9| . This connection is made 
concrete by considering the classical limit of a cavity pre- 
pared in a coherent state with « » 1. Neglecting the 
photon-number fluctuation, the single-electron spectrum then 

becomes Ej^(k) - ± -^(to - e^)^ + /2, with = yf{n)al 
being the classical Rabi frequency and (n) the average num- 
ber of photons. The bulk topology can be correspondingly 
characterized by Eq. ([3]), with 9k - arctan [Qt/((jj - 6^)1/2. 
We find P - \ for t < o) < t + and zero otherwise. 
Thus the topology of the electron system can be modified due 
to the interaction with the cavity. 

The preceding results were found within a single-particle 
picture, which in fact loses its precise meaning in the many- 
body Dicke-model regime we are actually dealing with. 
Namely, all the pseudospins simultaneously interact with the 
same cavity mode, and the exact many-body spectrum in- 
volves all of them. In order to simplify further progress, we 
notice that the interacting Hamiltonian admits a conserved 
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global quantity C - Z^iT^ + a^a. The simplest situation 
is when there are no photons in the cavity and C = -N/2 
reaches its lowest possible value. The combined eigenstate of 
the system is just the product state liA^t) - ® I0) where 
IU-) = li • ■ ■ i>, and the energy £0 = - ^k- For a finite 
number of photons in the cavity the problem becomes pro- 
gressively more involved, and reasonably simple results can 
only be found for n = 0,1,2, as well as for n » 1 (mean 
field), with n = C + N/2. 

We next wish to analyze the many-body ground state (at a 
fixed-n subspace) when the cavity is populated with photons. 
We focus only on the n - 1 case, as this already contains much 
of the relevant physics. In this case, the wave function can be 
expanded as 
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liAtot) - fllli) ® |1> + ^^*:IU-;Ta> ® |0) , (7) FIG.4. Main: The energy splitting of the two lowest-energy states, 



where the coefficients a and bk's are found by solving the 
Schrodinger equation "Kli^J^t) = Ei\t//l^^)\36\. The simplest 
case is when /i = 0, so that = f, and the ground-state energy 
is found to be £1 = -(A^ - l)t + a>/2 - ^jNY^ + (t - u)/2f, 
with Fi = 2i;(Q'j)^/A^. (For a finite ju 0, the result is more 
complicated, albeit similar in nature.) This means the ground- 
state energy is lowered at resonance compared to the non- 
interacting situation by the amount oc ^|N. If the number of 
pseudospins (or the strength of coupling Fi) exceeds a crit- 
ical value, the absolute ground state thus becomes populated 
with photons. For example, in the homogeneous coupling de- 
scribed above £1 = £0 for Fi - uP- /N at resonance. Al- 
though this transition is not innately topological, it can switch 
the topology of the system. Thus, at the critical coupling, the 
system switches to a different ground state populated by pho- 
tons, and the topology, which is related to the appearance of 
end modes, must be appropriately reexamined. To that end, 
however, we need to calculate the total many-body energy in 
a finite system, which is beyond the single-particle bulk rea- 
soning. 

Finite system. — The Majorana states are actually boundary 
modes, even though their existence can be inferred from the 
bulk spectrum. We use the discrete lattice model to numeri- 
cally diagonalize the Hamiltonian in the zero- and one-photon 
regimes, in order to explicitly identify zero modes. We have 
previously noted the existence of Majorana end modes of a 
free chain for ^ - Q. That result is extended to ju 0, af- 
ter finding an orthogonal transformation M such that Hid - 
MHioAf - 'ZLiWiri, withy^; = 2:^./;.r;: (here, 
j,m = 1, . . . , A^; p,r - 1,2), s^j(t,id) being the elements of 
the real orthogonal 2A^ x 2A^ matrix Ai [iTl. We can write also 
^^iD = 2m=i £m(c!„'c„, - 1/2), with ?;'„ = (j^ + i%)/2, and 
the eigenenergies e„ define the spectrum of the chain which 
contain the zero modes for \jj.\ < t. 

Focusing on the topological regime, we now analyze the 
evolution of the many-body zero mode as the first photon is 
starting to populate the cavity. We assume the cavity fre- 
quency to satisfy (f - |yu|) < to < 2{t - which allows 
for the resonant transitions only from the Majorana states to 



A£maj, for zero (C = 0) and one excitation (C = 1), as a function of 
the number of sites A'. For these plots, we used fi = 0.1, = 0.1, 
y6 = 0, and ai = 1 , all in units of /. Inset: The effective gap from the 
Majorana states to the continuum, Aiig^p, for C = 1, as a function of 
A' (using here a = 0.2 to exaggerate trends). 



the gapped fermionic continuum. Note that the parity of the 
system, p = y} •-.y^, is conserved, i.e., [p,'H] - 0, which 
means this is a good quantum number even in the presence of 
the cavity. The cavity may thus afifect the splitting between the 
parities, but not mix different parities. The interaction Hamil- 
tonian can be written in terms of the operators c^[,(ci„) instead 
of Til's which, in the rotating- wave approximation, becomes 
^int = I„n^mfAmf(cF,'cp)'cla + H.C., with A,„ the effective 
coupling and /(cp, c];.) being a linear function inF/7,F|,. As in 
the pseudospin model, here too we can use the number of ex- 
citations C = Yjn'^i^m + a' a as a conserved quantity for each 
parity p to find the spectrum and, as an example, we calculate 
the energy of the system for C = 1 . The corresponding eigen- 
values for the two parities are found from solving the eigen- 
value problem similarly to Eq. (|7]i, which boils down (see SM 
for details) to solving the equation 



2e„, 



0, 



(8) 



where e-t stand for the energy of the states of parity ±; ej^^ = 
Zje,n; r,r',p,p' = 1,2; and^-' = a/'-' Zjis'/^sf, - s';^s%) 
setting p - 0. The splitting between the lowest-energy 
states of different parities in the subspace C = 1 is denoted 
by A/imaj = - while the gap to the continuum by 
^£^gap = rnin (ef) - max (ef'*), with ^^^^'^^ being the ground- 
state (excited-state) energies. Note that such a definition can 
be generalized for any C > 1 , but not for C - 0. 

Figure |4] shows the energy splitting between the lowest en- 
ergy states, A/iniaj, for C = and C = 1, as well as the 
splitting (gap) between these two states and the continuum, 
A/igap for C - I. We see that in both cases the splitting be- 
tween the two states scales exponentially with the number of 
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sites, although slower for C = 1. Interestingly, the energy 
gap A/igap saturates at a value oc a for large A^. Physically, 
the resonant interaction brings the Majorana modes into the 
continuum spanned by the C = 1 states, but as this interac- 
tion increases, precisely one state per parity is pushed below 
the continuum and span a well-defined two-dimensional de- 
generate subspace. These degenerate modes are highly entan- 
gled electron-photon states, which we call Majorana polari- 
tons. Note that there is no Dicke-like enhancement y/Na of 
Aiigap for large (see Ei below Eq. (jTji), since it is exactly 
canceled by the overlap of the Majorana states with the contin- 
uum, which scales as 1 / yfN. Nevertheless, for large enough 
a, the Majorana polaritons are well defined. 

Summary. — To conclude, we analyzed the spectrum of a ID 
Kitaev chain coupled to a microwave cavity in both off- and 
on-resonant regimes. An off-resonant cavity can serve as a 
measuring device for the topological transitions, while in the 
resonant regime the topology can be changed depending on 
the photon number. We show that at the (mean-field) single- 
particle level the spectrum maps to that of classically-driven 
(Floquet) topological insulators. In the many-body picture, we 
find that the spectrum can exhibit a degenerate ground state 
even for a finite number of photons (despite strong hybridiza- 
tion of the original Majoranas with the continuum). These 
states are entangled electron-photon composites interpreted as 
Majorana polaritons. Our model can be relevant for Ising spin 
chains Il37l . one-dimensional nanowires with strong SOI in 
the presence of an i-wave superconductor (W\, as well as in 
the context of cold-atom physics |,19J . 
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CHAIN-PHOTON COUPLING IN SECOND-ORDER PERTURBATION THEORY 



In this Supplementary Material we give a detailed description of the second-order calculation for the two types of the wire- 
cavity coupling Hamiltonians. We will focus on the coupling via the chemical potential since the other type of coupling (via 
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the modification of f) turns out to be a straightforward extension. The BdG+cavity Hamiltonian in the k space reads 'H = 



^^BdG + ^^int + ^^ph: 



ph 



(9) 



where Uk - (-t sin k, 0, -it cos k + ji)) and - a + p cos A;, f > is the hopping parameter, // is the chemical potential, w is the 
intrinsic cavity frequency, and a (J3) are is wire-cavity coupling strengths via modification of the chemical potential (hopping 
parameter). We consider the off-resonant regime, oj » t+\fi\, and assume the frequency to be large enough for the interaction term 
to be treated perturb atively. The precise condition for the applicability of the perturbation theory will be discussed below. We 
use the Schrieffer- Wolff transformation to derive the second-order corrections. This implies finding an anti-Hermitian operator 
S such that 



'H^e^'He-^ = 'H + [S ,'H] + ^[S ,[S + . . . , 
is diagonal in the eigenbasis of H^dc + ^^ph- For the linear in a/i terms to vanish, we require = [^Bdc + ^^ph> '^l or 

where the (super)operator X. is defined formally by X.A = [H^dc + Hpi,,A], for an arbitrary operator A. This gives: 



S - /lim 



dsQk e 



.Sr ak\. -I- , 2fsin^ ... 

I > — iTr(fl'-fl) tAw + a) 



(10) 



(11) 



(12) 



disregarding terms of order u . This is seen by explicitly evaluating Tk{s)'- 



Tkis) = rfcos {2eks) + 



f sinA: 



sin(2eii)T-j. + 2- 



t cos k + fi (tsink ^ t cos k + ^ 



n sin (eks), 



where = In^^l. Inserting this ^(.(i) in the integral over time leads to the following result for S : 



k 



k k 



irUa' — a) — 



2fsinA: ,, 



rl(fl' + a) - i 



4(f cosfe + fi) 



rik ■ Tk(a' - a) 



ak_ 



.,. .2fsinA: . 
rAa' - a) + i rAa' + a) 



(13) 



(14) 



for w » f H- which allows us to neglect the l/aP terms. Now we are in position to derive the effective Hamiltonian within 
second order, Ai/int = [S, Hi„i]/2, which is given by: 



1 V-' ffifff L r 



-I- . + 2fsinA: ^, , j. 

(a -fl),T^(fl + a)] + i [Ti^Ao -t-fl),T^(a + a)] 



k,k 

akOk' _- . ^tal sin k 



Zj /,) l< k' /_i 



k,k' 



a'a + 



(15) 



This is the Hamiltonian used in the Main Text (MT). We identify two distinct contributions: the purely electronic first term 
that introduces many-body effects and the second term (dependent on the photon number) that renormalizes the p-wave pairing 
strength. Note that the off-diagonal terms oc a^[(a'')^+a^] were dropped here, since their subsequent Schrieffer- Wolff elimination 
would produce higher-order terms that are outside of our scope. 



MEAN-FIELD SOLUTION FOR FREQUENCY SHIFT 



In the MT, we argued that the second term in the Hamiltonian, Eq. (15 1, leads to a shift in the cavity frequency that can be 
measured, allowing to identify the transition point as well as the phase the ID Kitaev chain is in. The corresponding expression 
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for the frequency shift is 



4f X — 1 

6u-^y alsinkTl. (16) 



Performing then the substitution — > (t^), which neglects electronic fluctuations, we obtain a frequency shift {50)) (that is 
independent of the photon number if neglect the small backaction of this coupling on the electronic system). However, we did 
not explicitly evaluate the average (r^), and, moreover, did not discuss how the first term in Eq. ( 15 i was accounted for Note 
that the first term is oc Ijo), while the second oc l/w^, i.e., much smaller in the high-frequency regime of interest. The first term. 



ALT' V "l"'^!'' w 

^-^int = Zj ' (17) 



is also special in that it leads to interaction efifects. Let us first consider it separately. In a mean-field approximation, the 
many-body effects are eliminated by the substitution: 

w+^<4)-<^P<4)- (18) 

We assume, furthermore, that ai^ = a (i.e., there are no tunneling effects). This mean-field interaction can thus be accounted for 
by a chemical potential shift, namely yU — > jUgff - ^i- A with 

k 

where the average (. . .) is taken self-consistently over the renormalized states. Neglecting the frequency-shift term oc l/w^ in 
Eq. ( 15 1 (which could also be treated exactly if necessary), we can write the renormalized BdG Hamiltonian as follows: 

^BdG = - Z [(^ ^ + + ' sin krl] , (20) 

k 

whose eigenenergies are ef^ - ^fi + fj.^^ + Itji^f^ cos k. The expectation values for r'jf in the ground state of this Hamiltonian 
read: 

fsin^ tcosk + fi^g 

<^^> = ' <n> = — :s — ' (21) 



^k h 



which, in turn, enter in the self-consistent equation for A: 

.2 



= A > J- >A 

0) V 4 -J-n 



t COS k + u- A 
dk ^ . (22) 

^eiT 



Here, we transformed the sum into an integral, supposing a large number of sites, » 1 . In Fig. |5] we plot this A as a function 
of the chemical potential for several values of the dimensionless wire-cavity coupling parameter g = a~N/a)t. We see that 
A follows a linear dependence for t and it saturates for » f to values that are different depending on the coupling 

parameter g. Actually, we can obtain analytical results for A in the two limiting cases. For small <«; f, the asymptotic form is 
A - rjd with the slope: 

-- ^ - //eff = (1 - r)// = , (23) 



\+g iH-; 

which means the function A can be as large as // for large A^, but not exceed it. In this case, the cavity reduces the chemi- 
cal potential keeping the system in the topological phase, i.e., it stabilizes it by shifting the transition point. In the opposite 
(nontopological) case, » f, the function r becomes instead: 



For the perturbation theory to be internally consistent, we need to assume w » a VA^, which makes the second term in Eq. ( 15 1 
small compared to the hopping in //BdG- Concerning the first term in Eq. ( 15 i, it is of order (a y/N jojf'u <s; w if a VA'^ (w. 




FIG. 5. The cavity-induced chemical potential shift A as a function of the bare chemical potential ^l (in units of t) for different values of the 
effective coupling strength g = a-N jtoj. 



which can, nevertheless, be still comparable or even larger than t and/or A more stringent condition, g <K 1, would assure 
that Ai/int causes only small shifts in the electronic spectral properties. We can finally express the frequency shift Su) using the 
self-consistent wave functions to find: 



2^,2 



<(JW> 



sin^^ 



^fi + f?^ + 2fyUeff COS k 



2fig 



f 



dk- 



t sin^ k 



(25) 



f- + iTr. + Itj^efi COS k 



which reduces to the expression used in the MT. As can be seen from Fig. [5] for g <s: 1, we can approximate yUejf ^ fi. In the 
MT, we wrote {6a>} - 4Nt{a/cj)^r(fi/t) and plotted both r(x) and FaCx) = d^F/dx^. While the general expression for r(x) is 
too lengthy to be showed here, there are limiting cases when this becomes very simple. We will thus provide the expressions for 
r(ji:) and r2(x) for the regimes x <K 1, x ^ 1, and x » 1. For r(x), we obtain: 



r(x)- 



2 ' 1 ~ TT I ' for X <K 1 



37T 



(3-x), for 



1 



(26) 



1 

2x ' 



for X » 1 



which implies that {6(0)^^=0 - 2tN(a/oj) (maximum), and it scales as (dco) ~ l/i-i for » t. This means that {6aj) is just a 
fraction of the hopping parameter t. Note also that the function r(x) is continuos at x = 1, so there is no jump in the frequency 
shift. On the other hand, the limiting expressions for r2(x) read: 



r2(x) = 



1 



1 



1 + -X' 



for X « 1 



■ log |x - 1| , for X ^ 1 



(27) 



1 



for X » 1 



which diverges logarithmicaly as x 



ROTATING- WAVE APPROXIMATION FOR FINITE-CHAIN SPECTRUM 



In this section, we derive the rotating-wave approximation (RWA) Hamiltonian for a finite chain, as well as the resulting 
eigenenergy equation for C = 1 [Eq. (8) in the MT]. We first write the wire and interaction Hamiltonians in terms of Majorana 
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operators: 



N-l 



Hmt — i 



N-l N 



(a' + a) + H.c. 



(28) 



L M ./=i 
We can diagonalize the Hamiltonian in the absence of the cavity using the prescription in Ref. [1] of the MT and write 

N N 



m=l 



m=l 



(29) 



with e,„ > being the spectrum and 



( y\ 1 
y\ 




( y\ \ 


= M 


y\ 


% 




y'N 




^yl> 



(30) 



in terms of an orthogonal 2N x 2N matrix M (MM^ - M^M - I2n)- The fermionic operators are given by cl, - (y^, - iyf„)/2 
and Cm = (yl, + ;y^)/2. More specifically, we can write 



^ ' 

j=l p' = l,2 



(31) 



with the coefficients being the entries of the matrix M that can be found numerically. Note that if e„, is a solution, that -e 
is a solution too. Due to this symmetry, there are possible zero-energy solutions which, if existing, satisfy 

with the coefficients C''^ determined by the normalization conditions (not shown, see Ref. |[T|). Next, we write 



(32) 



m=l p=l,2 



(33) 



where Sjl^^ = r^^ j (since M is orthogonal). The spin-photon interaction Hamiltonian is then obtained as follows: 



/ 1 / J " / J 7+l,m j+l,m-^^ j,m' J,nr / j j,m j,m' 



r':,J'M +a) + H.c. 



= Z + (-l)""'c,„][?];, + (-l)"'-'c„,,](a^ +a) + H.c. , 



m,m' p,p' 



in in' 



N-l N 
" / J ' J+l,tn j+i,m'' j,m' j,m' / j j,m j,m' 



(34) 



Next we use the RWA to deal with the resonant regime thus excluding the fast-rotating terms. The coupling Hamiltonian in the 
interaction picture reads: 



(35) 



m,m' p,p' 
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When the frequency is such that it couples resonantly only the Majoranas to the continuum, but not the two bulk bands, the 
physics can be described in terms of a Majorana polariton, as outUned in the MT. Thus, we simplify the discussion by performing 
the RWA to the above Hamiltonian to get: 



H 



RWA 



m p,p' 

= Z Z « - <o ) (^0 + (4« + (-l/-'c„at) 

= ^ Am(cl,'c())cl,a + h.c. , 

m 

A„(?j.Fo) = Z « - <o ) (^0 + (-ir ^^o) . 



(36) 



Note that we labeled the zero modes by m = (for convenience and not meaning the first site). As mentioned in the MT, this 
Hamiltonian allows for the following conserved quantity: 



(37) 



which is easily checked to commute with the Hamiltonian. Next we compute the energy of the system for C — I (noting that 
C = is trivial). The ground state (GS) in the absence of the cavity is spanned by the following states: 



IC) = I0m)®|0^), 

l«/'f) = |lM)®|0^), 
iV-1 



(38) 



corresponding to the Majorana state being either empty (0) or occupied (1), while all the other gapped states are empty. The 
excited states are then written as 



|^i^) = |0M>®|0...1,-...0), 

|^i^) = |lM>®|0...1;...0>, 



(39) 



where j - 1,N - 1 . We are now in position to write down the general state corresponding in the presence of the cavity to C = 1 

(one excitation) as 



\^ex) = [(ao|OM> + aillM) ® 10 . . . 0)] ® \lph) + 



N-l 



Y,ibl\OM) + b{\lM)®\O..Aj...O) 

7=1 



® \Oph) . 



We need to choose the coefficients ao,i and ^ so that the state is an eigenstate of the total Hamiltonian: 



(Hm + H™"" + cua^a)|a>ex> = e|Oex> . 
First we apply to the state \%^y. 

^^SriOex) = Z ^" f^^ol^^^ + (-l)^'-'ai|OM» ® 10 . . . 1„ . . . 0)1 ® lOph) 

m,p,p' 

+ Z ^m'(-l)^"'[(^ollM) + (-l)^'-'^riOM))®|0...0)]®|lph). 
m,p,p' 

On the other hand, the effect of H\i, + //ph on this state reads: 



(40) 



(41) 



(42) 



(//lD+//ph)|Oex) 



cl> - } e, 



(aolOM> + oiIIm)) ® |0 . . . 0) ® \lph) 



e„il - 26j,Jibi\0M) + b{\lM))\0 ... 1; ... 0) 



® lOph) . 



(43) 
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2/i 
A 



9 



a 




hare Major anas 



dressed Majoranas 
(Majorana polaritons) 



FIG. 6. Sketch of the energy-level structure for the one-excitation regime (C = 1). The Majorana states in the absence of cavity (left) for parity 
-I- (red) and - (blue) are brought into the continuum via the interaction with the cavity. As a is increased, precisely one state for each parity 
is pushed below the continuum to form new Majoranas, now highly-entangled electron-photon states that can be called Majorana polaritons. 
Here, ji is the chemical potential, Aj stands for the initial gap from the zero mode to the continuum (in the absence of cavity), and w is the 
cavity frequency. 



We can now put everything together to find the equations for the energy: 

eao = (w - ^ e„,)flo + Yj i-^)'^"' ^'nf K , eb{ = a,, ^ sf - b{ ^ e„,(l - 26 j,^) , 

m m.p.p' p,p' m 

m tn,p,p' p,p' m 



(44) 



where it is evident that we are getting two decoupled sets of equations for each parity: (aq, {b"'}) and (oi , {fe™ )). One can easily 
obtain from these relations the eigenenergy equations: 



(1 - 25„,„,0 



= 0, 



(45) 



for the parity p - +{-), which are quoted in the MT. Besides the spectrum e we can also extract the set of coefficients 
(flo(i), {fe™Q))) (together with the normalization condition), which quantify the amount of photon-electron entanglement in the 
system. In the MT, we argued that the spectrum obtained by solving the Eq. (45 i yields one state for each parity emerging below 
the continuum of C = 1 states for sufficiently large coupling strenght a, which are degenerate (i.e., exponentially small splitting 
in the number of sites A^). These states, which we called Majorana polaritons, are shown schematically in Fig. [6] While for 
vanishingly small a the initially isolated Majoranas get mixed with the C - \ continuum of states, as a is increased further, new 
Majoranas reemerge, with well-defined parity and energy. 

I 



